
University of 
Hertfordshire 


THE NON-STELLAR INFRARED CONTINUUM OF 
SEYFERT GALAXIES 


A Alonso-Herrero 1 ’ 2 , A C Quillen 1 , C Simpson 3 , A Efstathiou 4 and M Ward 5 


1 Steward Observatory, Tucson 

2 University of Hertfordshire 

3 Subaru Telescope, Hilo 

4 Blackett Laboratory, Imperial College London 

5 University of Leicester 



Draft version November 29, 2000 

Preprint typeset using IAT^X style emulateapj v. 04/03/99 


THE NON-STELLAR INFRARED CONTINUUM OF SEYFERT GALAXIES a 

"THIS WORK IS BASED ON OBSERVATIONS COLLECTED AT UKIRT. THE UKIRT IS OPERATED BY THE 
JOINT ASTRONOMY CENTRE ON BEHALF OF THE UK PARTICLE PHYSICS AND ASTRONOMY RESEARCH 
COUNCIL. BASED ON OBSERVATIONS WITH THE NASA/ESA HUBBLE SPACE TELESCOPE, OBTAINED AT 
THE SPACE TELESCOPE SCIENCE INSTITUTE, WHICH IS OPERATED BY THE ASSOCIATION OF UNIVERSI- 
TIES FOR RESEARCH IN ASTRONOMY, INC. UNDER NASA CONTRACT NO. NAS5-26555. 

Almudena Alonso-Herrero 

Steward Observatory, The University of Arizona, Tucson, AZ 85721, USA 
Present address: University of Hertfordshire, Department of Physical Sciences, College Lane, Hatfield, Herts 

ALIO 9AB, UK 
Alice C. Quillen 

Steward Observatory, The University of Arizona, Tucson, AZ 85721, USA 

Chris Simpson 

Subaru Telescope, National Astronomical Observatory of Japan, 650 N. A‘Oh0k0 Place, Hilo, HI 96720, USA 

Andreas Efstathiou 

Astrophysics Group, Imperial College London, Blackett Laboratory, Prince Consort Rd, London SW7 2BZ, 

UK 

AND 

Martin J. Ward 

Department of Physics and Astronomy, University of Leicester, Leicester LEI 7RH, UK 
Draft version November 29, 2000 

ABSTRACT 

JHKL'M (1 — 5/xm) imaging of a sample of Seyfert 2 galaxies is presented. We have performed 
an accurate estimate of the near-infrared non-stellar nuclear fluxes. We confirm that the near-infrared 
nuclear continuum between 1 and 2.2 /xm of some Seyfert 2s is dominated by stellar emission, whereas 
the continuum emission at longer wavelengths (A = 3 — 5 /xm) is almost entirely non-stellar in origin. 
The non-stellar spectral energy distributions (SED) in the infrared (up to 15 /xm) of Seyfert galaxies 
show a variety of shapes, and they are well reproduced with the tapered disk models of Efstathiou & 
Rowan-Robinson (1995). We have used two models, one including an optically thin cone component 
found to fit the SED of NGC 1068, and a coneless model. Although our modelling of the SEDs does 
not allow us to favor either model to account for all the observed SEDs, we find that the viewing angle 
towards the central source is well constrained by both models. The galaxies in our sample have fitted 
values of the viewing angle in the range 0 V = 0° - 64°, for the assumed model parameters. We have also 
investigated non-stellar color-color diagrams (L' - M vs. H — M and V - M vs. H — V). The colors 
of the Seyfert galaxies with viewing angles 6 y < 30° are better reproduced with the cone model. These 
diagrams provide a good means to separate Seyfert 2s with moderate obscuration ( Ay <; 20 mag from 
hard X-ray observations) from those with high obscuration. 

The ground-based 4.8 /xm and ISO 9.6 /xm luminosities are well correlated with the hard X-ray lu- 
minosities of Seyfert Is and 2s. These continuum emissions appear as a good indicator of the AGN 
luminosity, at least in the cases of hard X-ray Compton-thin Seyfert galaxies (Nr < 10 24 cm~ 1 2 ). We 
finally stress the finding that some Compton thick galaxies show bright non-stellar emission at 5 /xm. 
This suggests that the near-infrared emission in Seyfert galaxies is produced in an extended component 
illuminated by the central source, that is more visible from all viewing angles, providing a good expla- 
nation for the differing Nh/Av ratios found in some Seyfert 2s. We discuss possible implications of 
mid-infrared surveys for the search of counterparts of highly obscured hard X-ray sources. 

Subject headings: Galaxies: nuclei - Galaxies: Seyfert - Galaxies: active galaxies: photometry - 
galaxies: stellar content - infrared: galaxies 


1. INTRODUCTION 

There is now considerable evidence that the infrared 
emission in active galaxies (Seyfert galaxies and radio quiet 
quasars) is thermal in origin, namely hot dust emission 
(e.g., Rieke & Low 1975; Rieke 1978; Barvainis 1987; 
McAlary & Rieke 1988; Sanders et al. 1989; Danese et 
al. 1992 and recently Andreani, Franceschini, & Granato 
1999; Haas et al. 2000; Polletta et al. 2000 and refer- 


ences therein). The strong mid and far infrared emission 
in active galactic nuclei is often interpreted as produced 
by a disk-like (Sanders et al. 1989) or torus-like (e.g., 
Pier & Krolik 1992; Granato &; Danese 1994; Efstathiou 
& Rowan-Robinson 1995 among others) dust configura- 
tion heated by the central source. This obscuring material 
which hides a direct view of the central source in type 2 
objects, absorbs a large fraction of the X-ray /UV/«pticaI 
energy emitted by the central source and reradiates it in 
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the infrared spectral range. 

The nature of the non-stellar infrared continuum in 
Seyfert 2 galaxies where the stellar component dominates 
the nuclear near-infrared emission (Edelson & Malkan 
1986; Alonso-Herrero, Ward, & Kotilainen 1996), has 
eluded us due to the difficulty in estimating the non-stellar 
contribution at near-infrared wavelengths (up to 2/xm). 
High resolution near-infrared images observed with NIC- 
MOS on the HST have now been used to detect and mea- 
sure the nuclear non-stellar emission in Seyfert galaxies 
(Quillen et al. 2000b). This work reported that almost 
100% of all Seyfert 1-1.9 galaxies and 50% of all Seyfert 
2s show point sources at 1.6 /im. This unresolved emis- 
sion is variable in some of the Seyfert 2s and thus is pre- 
sumably non-stellar emission associated with the central 
engine (Quillen et al. 2000a). The high detection rate 
of point sources in the near-infrared is in contrast with 
the small number of point sources reported in Seyfert 2s 
in the HST optical survey of Malkan, Gorjian, & Tam 
(1998). Previous attempts to model the non-stellar spec- 
tral energy distributions (SEDs) of Seyfert 2s (e.g., Edel- 
son & Malkan 1986, and recently Fadda et al. 1998, and 
references therein) suffered from a number of problems, 
including large aperture photometry, and uncertain deter- 
minations of the stellar component which may dominate 
the near-infrared emission up to 2.2 /im of Seyfert 2s. To 
date, a detailed modelling of SEDs has been restricted to 
individual galaxies (e.g., NGC 1068 Efstathiou, Hough, &; 
Young 1995; NGC 3281 Simpson 1998; CenA Alexander 
et al. 1999; the Circinus galaxy Ruiz et al. 2000). An 
additional problem may exist when non-simultaneous ob- 
servations are used since some Seyfert 1.8 and 1.9 s are now 
found to show variability at near-infrared wavelengths over 
time scales of several months (Quillen et al. 2000a from 
J/ST/NICMOS observations at 1.6 /zm). 

Models of dusty tori with different geometries (e.g., Pier 
& Krolik 1992; Granato & Danese 1994; Efstathiou & 
Rowan-Robinson 1995; Granato, Danese, & Ftanceschini 
1997) predict bright infrared emission, making the L-band 
(at 3.5 /zm) and M-band (at 4.8 /zm) well suited to study 
the nature of the non-stellar emission in Seyfert galaxies. 
Because of the strong dependence of the infrared emission 
with the viewing angle and the geometry of the obscur- 
ing material, infrared SEDs can be used to constrain the 
torus parameters in Seyfert galaxies. An additional ad- 
vantage is that the stellar contribution is greatly reduced 
longwards of 2 /zm, even for Seyfert 2s (Alonso-Herrero et 
al. 1996). The use of observations at longer wavelengths 
> 50 /im (e.g., IRAS and ISO measurements), whilst useful 
in determining the bolometric emission, are problematic 
when trying to separate the contributions from the cen- 
tral AGN and a possible circumnuclear star burst, because 
of the large apertures employed. 

In this paper we present near-infrared (from 1 to 5 /zm) 
imaging observations of a small sample of Seyfert 1.9- 
2 galaxies which are complemented with ifST/NICMOS 
and ISO imaging, covering a spectral range from 1 up to 
15 /xm. In a future paper (Quillen et al. 2001, in prepara- 
tion) we will discuss the infrared properties and SEDs of 
the Seyfert 1.8-2s in the CfA sample. The paper is orga- 
nized as follows, in Section 2 we describe the observations, 
in Section 3 we determine the non-stellar emission, and in 
Section 4 we analyze the observed non-stellar energy dis- 


tributions, non-stellar color-color diagrams, and compare 
them with the outputs from torus models. In Section 5 
we will show that the 5 — 10 /zm fluxes are a good indi- 
cator for the AGN luminosity for Compton-thin AG Ns. 
Conclusions are presented in Section 6. 

2. OBSERVATIONS 
2.1. UKIRT Observations 

We obtained J (A c = 1.25 /zm), H (A c = 1.65 /zm), K 
(Ac = 2.20 /zm), V (Ac = 3.80 /im) and M (A c = 4.80 /zm) 
imaging of a sample of fourteen Seyfert 1.9 and Seyfert 2 
galaxies (see Table 1) with the infrared camera IRCAM3 
on the 3.9-m United Kingdom Infrared Telescope (UKIRT) 
during two observing runs in September 1997 and April 
1998. The size of the J, H and K images is 256 x 256 pixels 
with pixel size of 0?143 pixel -1 , whereas the size of the 
V M images is 64 x 64 pixels, with pixel size 0^/286 pixel - 1 . 

Standard reduction procedures were applied. Condi- 
tions were photometric during both runs, so standard star 
observations from the Elias et al. (1982) list (for JHKL f ) 
and from the UKIRT list of photometric standards (for M) 
were used to perform the photometric calibration of the 
images. The major source of uncertainty in near-infrared 
observations comes from the background subtraction, es- 
pecially at the longer wavelengths (1/ and M). Typical er- 
rors from the photometric calibration are 0.09, 0.04, 0.08, 
0.07 and 0.10 mag in J, H , K , V and M, respectively. 
The full width half maximum (FWHM) seeing was mea- 
sured from the standard stars 0.6 — 0.7" at K band for 
both runs. 

Standard synthetic aperture photometry was performed 
on all the images using a ^-diameter circular aperture 
(Table 2). The errors associated with the sky subtraction 
at L and M are given together with the aperture photom- 
etry in the last two columns of Table 2. 

2.2. HST/NICMOS Observations 

We searched the HST archive for broad-band NIC- 
MOS observations of our sample of galaxies. In Table 1 
columns (4) and (5) we list the camera and filter, and the 
observation date respectively. In columns (6) and (7) we 
indicate if variability has been found, and the time scale 
from Quillen et al. (2000a). The plate scales for cam- 
eras NIC1 and NIC2 are 0.045” pixel -1 and 0.076” pixel -1 
respectively. 

The images were reduced with routines from the pack- 
age NicRed (McLeod 1997). The main steps in the data 
reduction involve subtraction of the first readout, dark cur- 
rent subtraction on a readout-by-readout basis, correction 
for linearity and cosmic ray rejection (using fullfit), and 
flatfielding. Darks with sample sequences and exposure 
times corresponding to those of our observations were ob- 
tained from other programs close in time to ours. Usually 
between 10 and 20 darks were averaged together (after 
the subtraction of the first readout) for a given sample se- 
quence. Flatfield images were constructed from on-orbit 
data. 

The photometric calibration of the NICMOS images was 
performed using the conversion factors based on measure- 
ments of the standard star P330-E (Marcia Rieke 1999, pri- 
vate communication). Comparisons of the ground-based 
if-band aperture photometry with the NICMOS F160W 
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Table 1 

Log of the UKIRT and HST/NICMOS Observations. 


Galaxy 

(i) 

UKIRT 

(2) 

Date 

(3) 

NICMOS 

U) 

Date 

(5) 

var 

(6) 

timescale 

(7) 

NGC 1052 

JHKL'M 

Sep 1997 

NIC2 F160W 

Sep 1998 



NGC 1068 

JHKL'M 

Sep 1997 

NIC2 FI 10W, F160W, F222M 

Feb 1998 

no 

??? 

NGC 1097 

JHKL'M 

Sep 1997 



NGC 2992 

JHKL'M 

Apr 1998 

NIC2 F205W 

Nov 1998 



NGC 4968 

JHKL'M 

Apr 1998 

NIC2 F160W 

May 1998 



NGC 5252 a 

M 

Apr 1998 

NIC1 FI 10W, F160W 

Mar 1998, Apr 1998 

no 

1 month 

NGC 5506 
NGC 7172 

JHKL'M 

JHKL'M 

Apr 1998 
Sep 1997 

NIC2 F160W, F205W 

Apr 1998, Aug 1998 

NGC 7217 

JHKL'M 

Sep 1997 

NIC1 F110W, F160W 

Aug 1997 



Mkn 1 

JHKL' 

Sep 1997 

NIC1 F160W 

Sep 1997 



Mkn 348 ft 

M 

Sep 1997 




Mkn 533 

JHKL'M 

Sep 1997 

NIC1 FI 10W, F160W 

Sep 1998, Sep 1997 

yes 

14 months 

Mkn 573 a 

M 

Sep 1997 

NIC1 F110W, F160W 

Aug 1998 

no 

12 months 

MCG -5-23-16 

JHKL'M 

Apr 1998 


Note. — a JHKL' imaging can be found in Alonso-Herrero et al. (1998). 


Table 2 

UKIRT PHOTOMETRY THROUGH A 3"-DIAMETER APERTURE FOR THE SEYFERT 2 SAMPLE. 


Galaxy J H 


NGC 1052 

11.40 

10.66 

NGC 1068 

10.42 

9.08 

NGC 1097 

11.95 

10.99 

NGC 2992 

12.34 

11.50 

NGC 4968 

13.26 

12.39 

NGC 5252 a 

13.27 

12.56 

NGC 5506 

12.09 

10.46 

NGC 7172 

12.74 

11.76 

NGC 7217 

11.99 

11.27 

Mkn 1 

14.12 

13.33 

Mkn 348 a 

13.52 

12.72 

Mkn 533 

13.38 

12.28 

Mkn 573 a 

12.87 

12.27 

MCG -5-23-16 

11.94 

11.01 


K 

V 

M 

10.39 

9.77 ±0.02 

9.08 ±0.23 

7.54 

4.53 ± 0.02 

3.24 ± 0.02 

10.75 

10.19 ±0.17 

> 10.50 

10.94 

10.01 ± 0.10 

9.15 ±0.30 

11.61 

9.98 ± 0.05 

8.65 ± 0.20 

12.03 

10.58 ± 0.04 

10.04 ±0.30 

9.14 

7.11 ± 0.03 

6.22 ± 0.02 

10.91 

9.50 ± 0.02 

8.56 ± 0.05 

10.96 

10.88 ±0.15 

> 10.60 

12.89 

11.26 ±0.10 


12.17 

10.50 ± 0.04 

9.06 ±0.10 

11.21 

9.14 ± 0.02 

7.95 ± 0.08 

11.78 

10.14 ±0.04 

8.98 ± 0.30 

10.22 

8.51 ±0.02 

7.67 ±0.06 


Note. a JHKL '- band photometry from Alonso-Herrero et al. (1998). 

The quoted errors for the V and M-band photometry account only for the sky subtraction uncertainties, except 
for the Alonso-Herrero et al. (1998) L'-band photometry which are the photometric errors. 

The Af-band limits for NGC 1097 and NGC 7217 are 3 a. 
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Table 3 

Mid-infrared ISO fluxes. 


Galaxy 

(i) 

A 

(2) 

flux 

[mJy] 

(3) 

reference 

(4) 

NGC 1068 

5.9 /x m 

8660 

Rigopoulou et al. (1999) 


7.7 /xm 

14500 

Rigopoulou et al. (1999) 

NGC 3227 

6.75 /xm 

294, 249 

Clavel et al. (2000) 


9.63 /xm 

382, 372 

Clavel et al. (2000) 


16 /xm 

1220 

Perez Garcia & Rodriguez Espinosa (2000) 

NGC 3281 

11.4 /xm 

530 

this work 

NGC 4151 

5.9 /xm 

882 

Rigopoulou et al. (1999) 


7.7 /xm 

1090 

Rigopoulou et al. (1999) 


16 /xm 

4120 

Perez Garcia & Rodriguez Espinosa (2000) 

NGC 5252 

6.75 /xm 

14.8 

this work 


9.63 /xm 

27.1 

this work 


11.4/xm 

34.6 

this work 


15 /xm 

44.8 

this work 

NGC 5506 

5.9 /xm 

646 

Rigopoulou et al. (1999) 


7.7 /xm 

713 

Rigopoulou et al. (1999) 


11.4 /xm 

890 

this work 

NGC 5548 

6.75 /xm 

170 

Clavel et al. (2000) 


9.63 /xm 

270 

Clavel et al. (2000) 


16 /xm 

440 

Perez Garcia & Rodriguez Espinosa (2000) 

NGC 7172 

6.75 /xm 

160 

this work 


15 /xm 

240 

this work 

NGC 7469 

5.9 /xm 

270 

Rigopoulou et al. (1999) 


6.75 /xm 

433 

this work 


7.7 /xm 

539 

Rigopoulou et al. (1999) 


15 /xm 

1290 

this work 

Mkn 1 

5.9 /xm 

20.9 

Rigopoulou et al. (1999) 


7.7 /xm 

88.3 

Rigopoulou et al. (1999) 

Mkn 348 

6.75 /xm 

65 

this work 


15 /xm 

270 

this work 

Mkn 533 

6.75 /xm 

259, 213 

Clavel et al. (2000) 


9.63 /xm 

345 

Clavel et al. (2000) 

IC 4329A 

6.75 /xm 

591, 477 

Clavel et al. (2000) 


9.63 /xm 

890, 739 

Clavel et al. (2000) 
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filter show that differences are always less than 15%, 
whereas the comparison of ground-based #-band photom- 
etry with that of the NICMOS F222M and F205W filters 
agrees within 25%. 

2.3. Comparison sample of Seyfert Is 

As a comparison sample, we also searched the HST 
archive for Seyfert 1-1.5 galaxies with 1 — 2.2^m NIC- 
MOS observations as well as with ground-based L (or V ) 
and M small aperture photometry and ISO observations. 

2.4. ISO Observations 

We searched the literature for available mid-infrared 
ISO (from 5 pm up to 15 pm) continuum fluxes. In addi- 
tion, for those galaxies with unpublished data, we obtained 
images from the ISO archive. The photometric calibration 
was that of the pipeline for which the typical errors are 
±15%. The ISO fluxes were measured using 12"-diameter 
apertures and aperture corrections measured from model 
PSFs. This procedure is very similar to that employed 
by Clavel et al (2000). NGC 5506 and NGC 7469 looked 
like point sources so contamination by nearby star forma- 
tion was not likely to be a large problem. For the bright- 
est source (NGC 5506) some additional problems may be 
caused by saturating the detector. The wavelengths, fluxes 
and references for the ISO data for the galaxies in our sam- 
ple of Seyfert 2s and Seyfert Is are given in Table 3. 

3. DETERMINATION OF THE NON-STELLAR FLUXES 
3.1. The 1 — 2.2 pm spectral range 

As mentioned in the introduction it is essential to ob- 
tain very accurate estimates of the non-stellar emission to 
understand and model the SEDs of Seyfert 2s. The stel- 
lar component may dominate most of the nuclear flux for 
wavelengths up to 2 pm in Seyfert 2s (Alonso-Herrero et 
al. 1996). The relatively small field of view of the ground- 
based JHK (36" x 36") and NICMOS images (19" x 19" 
and 11" x 11" for NIC2 and NIC1 observations respec- 
tively) prompted us to try a nucleus + bulge decomposi- 
tion of the observed radial surface brightness profiles. We 
used a similar method to that employed in Alonso-Herrero 
et al. (1996). For the ground-based UKIRT observations 
the nuclear component was assumed to be a delta func- 
tion convolved a Gaussian profile to account for the see- 
ing effects (see Alonso-Herrero et al. 1996 for more de- 
tails). For the NICMOS observations, the nuclear compo- 
nent was represented with a surface brightness profile of a 
point spread function (PSF) generated with the Tiny Tim 
software (Krist et al. 1998). The bulge component was 
represented with a de Vaucouleurs r* profile. The surface 
brightness profiles were extracted out to radial distances 
of 5.7", 4.6", and 2.6" for the ground-based data, NIC2 
and NIC1 data respectively. 

In general, we performed a nucleus ± bulge deconvo- 
lution to the radial surface brightness profiles, as advo- 
cated by Simpson (1994). For Mkn 348, however, only 
lower angular-resolution data were available (from Alonso- 
Herrero et al. 1998), and we opted for the method of Simp- 
son (1998), which involves scaling and subtracting the J- 
band image to the H and K images, according to the 
counts in a circumnuclear annulus. This method assumes 
that the non-stellar contribution at J is negligible and, 


like Simpson we confirmed this by extrapolating 

the stellar HkL*m fluxes to the J-band. 

After performing the nucleus ± bulge deconvolution we 
measured the unresolved flux, which will be assumed to be 
non-stellar in origin, and its contribution within a 3-arcsec 
diameter aperture. The non-stellar fractions within this 
aperture are given in columns (7), (8) and (9) for JHK 
respectively for the ground-based data, and in parenthe- 
sis for the NICMOS data when available. The non-stellar 
fraction within a 3-arcsec aperture (which corresponds to 
projected angular sizes of between 200 pc and 1.5 kpc for 
our sample of galaxies) gives an estimate of the domi- 
nance or otherwise of the non-stellar component at a given 
wavelength. In Table 4, columns (2), (3) and (4) we give 
the JHK non-stellar fluxes using the NICMOS non-stellar 
fractions (when available) and the ground-based photome- 
try. Note that we did not detect an unresolved component 
at 1 - 2.2 pm in NGC 1097 and NGC 7217, and therefore 
these two galaxies are excluded from the following analysis. 
We performed a similar analysis on the NICMOS images 
as with the Seyfert 2s to derive the non-stellar fluxes from 
1 to 2.2 pm (see Table 5). 

As discussed in Alonso-Herrero et al. (1996) the main 
source of uncertainty in the determination of the unre- 
solved fluxes from the ground-based data is the fitted value 
to the seeing. Because of the relatively small field of view 
of the images, no stars are present in the images. How- 
ever, the fitted values of the seeing from the profile de- 
convolution are very similar to the values of the FWHM 
measured from the standard stars observed close in time. 
On the other hand, the advantage of using NICMOS im- 
ages is that the NIC1 and NIC2 PSF is very stable, is very 
well sampled and can be well modeled with the TinyTim 
software. Therefore the estimated non-stellar fluxes from 
NICMOS profiles are less affected by uncertainties. An 
upper limit to the errors associated with the non-stellar 
flux determination can be obtained by comparing the non- 
stellar fractions derived from the ground-based and the 
NICMOS images. The differences are never greater than 
45%, and tend to increase for shorter wavelengths, where 
the non-stellar contributions are smaller. 

Our estimates are in relatively good agreement (to 
within the uncertainties given above) for those galaxies in 
common with Zitelli et al. (1993), Kotilainen et al. (1992), 
Alonso-Herrero et al. (1996) and Quillen et al. (2000b). 
We confirm the findings of Alonso-Herrero et al. (1996), 
that the 1 to 2.2 ^m continuum of some Seyfert 2 galax- 
ies is dominated by stellar emission. However, the stellar 
contribution decreases significantly at longer wavelengths. 

3.2. The 3 — 15 pm spectral range 

At the L'-band and especially at the M- band the host 
galaxy does not make a strong contribution. We used the 

V non-stellar fractions derived by Alonso-Herrero et al. 
(1998) for the galaxies in common with our study. For 
the rest of the sample we estimated the non-stellar frac- 
tions at V scaling the peak of the PSF surface brightness 
profile to that of the galaxy, and assumed that the PSF 
integrated flux represents the unresolved component (see 
Alonso-Herrero et al. 1998 for more details). An alter- 
native way of estimating the stellar contribution in the 

V can be done by using the typical colors K - L! of a 
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Table 4 

Non-stellar fluxes and non stellar contribution within a 3"-diameter aperture for the 

Seyfert 2 sample. 


Galaxy 

(i) 

f(J) 

[mJy] 

(2) 

/( H) 
[mJy] 
(3) 

f(K) 

[mJy] 

(4) 

f(L') 

[mJy] 

(5) 

/(M) 

[mJy] 

(6) 

J 

[%} 

(7) 

H 

[%] 

(8) 

K 

[%] 

(9) 

V 

[%] 

(10) 

NGC 1052 

no 

no 

<0.9 

<20.2 

38.0 

no 

no (0) 

< 2 

< 65 

NGC 1068 

9.8 

97.6 

449.6 

3690.8 

8245.0 

17 (9) 

47 (41) 

69 (71) 

95 

NGC 2992 

no 

< 1 

2.8 

22.7 

35.7 

no 

no 

10 (23) 

91 

NGC 4968 

no 

0.6 

3.7 

23.2 

56.5 

no 

5(7) 

25 

90 

NGC 5252 

no 

0.7 

1.0 

11.1 

15.7 

no 

4(7) 

10 

75 

NGC 5506 

13.8 

59.0 

120.4 

340.1 

530.0 

59 

79 (89) 

64 (83) 

95 

NGC 7172 

no 

< 0.4 

3.4 

30.0 

61.4 

no 

<2 

12 

75 

Mkn 1 

< 0.07 

0.15 

0.8 

6.4 

... 

<2 

3(7) 

18 

80 

Mkn 348 

< 1 

0.6 

2.7 

14.5 

38.7 

<2 

7 

30 

91 

Mkn 533 

1.0 

5.0 

12.3 

53.0 

108.6 

23 (14) 

31 (40) 

57 

95 

Mkn 573 

0.2 

0.6 

3.2 

18.8 

41.3 

(2) 

(6) 

25 

85 

MCG -5-23-16 

1.1 

3.7 

10.7 

79.5 

139.4 

4 

9 

20 

80 


Note. — The zero points used for JHKL'M are: 1600 Jy, 1020 Jy, 657 Jy, 252 Jy and 163 Jy respectively. The 
3"-diameter M-band fluxes are assumed to be all non-stellar. “No” means that no unresolved component was 
detected. 


Table 5 

Non-stellar fluxes and non stellar contribution within a 3"-diameter aperture for the 

Seyfert 1 sample. 


Galaxy 

a) 

/( F110W) 
[mJy] 
(2) 

/(F160W) 

[mJy] 

(3) 

/(F222M) 

[mJy] 

(4) 

m 

[mJy] 

(5) 

f(M) 

[mJy] 

(6) 

J 

[%} 

(7) 

H 

[%] 

(8) 

K 

[%] 

(9) 

NGC 3227 

• * • 

10.6 

22.6 

78.3 ± 4.4 

72 ±27 

... 

21 

32 

NGC 4151 

69.0 

103.6 

177.5 

325 ± 10 

449 ± 34 

82 

81 

77 

NGC 5548 

. . . 

15.0 

31.6 

98.6 ± 3.7 

100 ± 21 

. . . 

57 

59 

NGC 7469 

16.2 

39.0 

67.8 

159 ±5 

259 ± 33 

45 

57 

63 

IC 4329A 

... 

49.9 

101.8 

210 ±5 

170 ± 20 

‘ • • 

71 

75 


Note. — The L and M fluxes are from Ward et al. (1987), and are assumed to be all non-stellar in origin. 
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normal galaxy (Willner et al. 1984, and assuming a color 
transformation L — V = 0.2, Ward et al. 1982) and the 
stellar fractions at the ff-band. We find that the stellar 
contribution within a 3"-diameter aperture at the V are 
always less than 25% (except NGC 1052 for which we find 
a 65% of the total 3"emission is stellar emission) and are in 
good agreement with the values estimated from the imag- 
ing data. The same procedure was applied to the M - band 
fluxes using the K - M color of a normal galaxy (Rieke & 
Lebofsky 1978), and estimated that the stellar contribu- 
tion at 4.8 /xm within a 3"-diameter aperture is always less 
than 10%. This is supported by the compact appearance 
of all the observed sources at this wavelength. In Table 4, 
columns (5) and (6) we list the non-stellar fluxes at V and 
M. 

The small aperture L and M fluxes for the Seyfert 1 
sample (compiled by Ward et al. 1987) are assumed to 
be all non-stellar in origin. The non-stellar fluxes along 
with the non-stellar fractions within a 3-arcsec diame- 
ter aperture are given in Table 5. Note that the non- 
stellar 1.2 — 2.2 /xm fluxes are not simultaneous with the 
3.S-4.8 \i m observations, which may add some uncertainty 
to the observed SEDs. 

4. THE NON-STELLAR INFRARED EMISSION OF SEYFERT 
GALAXIES 

Because the non-stellar contribution to the near-infrared 
continuum emission from Seyfert Is usually dominates over 
the stellar emission, unlike the situation for Seyfert 2s (see 
e.g., Edelson & Malkan 1986; Kotilainen et al. 1992 and 
Alonso-Herrero et al. 1996), the SEDs of the non-stellar 
component is better defined in the Seyfert Is. It is agreed 
that the infrared (up to 100 /xm) SEDs of Seyfert Is are well 
represented with a power law (S u oc i/~ a ) with a canonical 
spectral index a ~ 1 (e.g., Edelson & Malkan 1986; Ward 
et al. 1987; Edelson, Malkan, & Rieke 1987; Granato & 
Danese 1994; Clavel et al. 2000). Similar spectral indices 
have been derived from the infrared SEDs of quasars (Edel- 
son 1986; Neugebauer et al. 1987; Polletta et al. 2000; 
Haas et al. 2000). Fadda et al. (1998) have compiled 
infrared data from the literature for a sample of Seyfert 
Is and Seyfert 2s, and found that the non-stellar SEDs of 
Seyfert 2s are always steeper than those of Seyfert Is, in 
line with the finding that the 1 — 2.2 /xm non-stellar fluxes 
of Seyfert 2s are significantly smaller than those of Seyfert 
Is. However, the work of Fadda et al. (1998) and that of 
Murayama, Mouri, & Taniguchi (2000) suffer from a num- 
ber of limitations when attempting to place constrains on 
the torus geometry based on the non-stellar infrared SEDs 
or colors. These works make use of non-simultaneous ob- 
servations, observations from different sources, uncertain 
non-stellar fluxes, and relatively large aperture N - and Q- 
band or IRAS fluxes which may be contaminated by other 
contributions. Moreover, in order to understand the na- 
ture of the non-stellar SEDs of Seyfert 2 galaxies and in 
turn to constrain the torus geometry and derive the in- 
frared extinctions to the central regions and their relation 
with the hard X-ray absorption columns, one cannot av- 
erage the non-stellar SEDs of Seyfert 2s (or Seyfert Is). 

An additional problem often present in past statistical 
studies of SEDs is the grouping of galaxies into the Seyfert 
1 and Seyfert 2 categories. A galaxy may be classified as 
a Seyfert 2 when it is observed in the optical (that is, no 


broad lines), but it may still show broad lines in the in- 
frared, where the extinction is greatly reduced (see for in- 
stance, Ruiz, Rieke, &; Schmidt 1994; Goodrich, Veilleux, 
& Hill 1994; Veilleux, Goodrich, & Hill 1997). If using the 
latter observations, one would classify it as a Seyfert 1.8 or 
Seyfert 1.9. A better approach to the understanding of the 
infrared SED of Seyfert galaxies is to study each galaxy 
separately and determine whether a single torus geometry 
is able to fit most of the galaxies. Although in our study 
we do not claim a complete or unbiased sample of Seyfert 
Is and 2s, because we have estimated accurate non-stellar 
SEDs, we can compare these with current torus models 
and see whether there is a general agreement across the 
range of properties present in this sample. We will make 
use of Efstathiou & Rowan-Robinson (1995) models, and 
discuss other models available in the literature. In addi- 
tion to the galaxies presented in this work, we will fit the 
non-stellar infrared SED of NGC 3281 (non-stellar fluxes 
from Simpson 1998). 

4.1. The Efstathiou & Rowan-Robinson (1995) models 

For our analysis we consider the tapered disk model of 
Efstathiou & Rowan-Robinson (1995). In this geometry, 
the height of the disk increases with the radial distance, 
but tapers off to a constant height in the outer part. We 
will use a simple torus model (coneless) and a compos- 
ite model with a torus and an optically thin cone. The 
latter model gave a better fit to the near-infrared emis- 
sion in NGC 1068 (Efstathiou et al. 1995). The reasons 
for adding the optically-thin component rather than using 
a simple torus model to fit the SED of NGC 1068 were 
twofold. First, this additional optically thin cone com- 
ponent increases the continuum around 3 — 5/xm, a con- 
dition necessary to fit the almost flat (in v S v ) observed 
SEDs of Seyfert Is (see also below). This problem seems 
to affect all the current torus models (e.g., Pier & Krolik 
1993; Granato & Danese 1994, and the model used here). 
Second, the simple torus model predicts strong absorption 
features at 9.7 /xm for Seyfert 2s, which are not usually 
observed (Roche et al. 1991). See also the discussion in 
Efstathiou et al. (1995). 

The model assumes a maximum dust temperature (dust 
sublimation temperature) of T = 900 K, and the following 
parameters: ri/r 2 = 0.01 (ratio of the inner and outer 
radii), h/r 2 =0.1 (ratio of the height scale and the outer 
radius) and equatorial UV optical depth of ruv = 1200, 
with a half-opening angle of the toroidal cone of 6 C = 30°. 
We will also explore for comparison purposes another 
coneless model with a dust sublimation temperature of 
T = 1500 K. 

A beaming factor was introduced by Efstathiou et al. 
(1995) to basically to boost the near-infrared flux, as the 
inclusion of conical dust by itself did not provide sufficient 
emission in this spectral range. The beaming factor is 
defined as the ratio of the radiation intensity emitted by 
the continuum source towards the optically thin cone with 
respect to that directed towards the torus, 

t _ Sy{0 £ &c) 

/b ~ s v (0 >e c ) K J 

where 9 C is the semi-opening angle of the torus, and 6 is 
the angle measured from the pole. The beaming factor for 
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Table 6 

Viewing angles, UV optical depths, and 

EQUIVALENT VISUAL EXTINCTIONS FOR THE 
MODELS. 


<?v 

n 

njv 

(mag) 

A v 

(mag) 

90 

1200 

221 

84 

1197 

221 

79 

1018 

188 

73 

914 

169 

64 

812 

150 

50 

700 

129 

44 

658 

121 

40 

620 

114 

36 

558 

103 

30 

263 

49 

30 

50 (24) 

9(4) 

25 

2.5 (0) 

0.5 (0) 

20 

2.5 (0) 

0.5 (0) 

10 

2.5 (0) 

0.5 (0) 

0 

2.5 (0) 

0.5 (0) 


Note. — In columns (2) and (3) in brackets are the 
values for the T — 1500 K model. 


Table 7 

Infrared non-stellar SED modelling. 


Galaxy 

n 

X 2 /" 

e v 

model 

NGC 1068 

7 

1.0- 1.6 

40° 

Cone, T = 900 K and No cone, T = 1500 K 

NGC 3227 

7 

1.0 

30° 

No cone, T = 900 K 

NGC 3281 

6 

CN 

CS 

1 

CO 

H 

40° 

Cone and No cone, T = 900 K 

NGC 4151 

8 

0.8 

10° 

Cone, T = 900 K 

NGC 5252 

8 

1.6 

30° 

Cone, T = 900 K 

NGC 5506 

7 

1.3 

30° 

No cone, T = 1500 K 

NGC 5548 

7 

0.9 - 1.0 

10° - 30° 

Cone, T = 900 K and No cone, T = 1500 K 

NGC 7172 

6 

0.7 -1.6 

40° 

Cone, T = 900 K and No cone, T = 1500 K 

NGC 7469 

9 

1.1 

o 

O 

CO 

1 

o 

O 

Cone, T = 900 K and No cone, T — 1500 K 

Mkn 1 

6 

1.2 

40° - 64° 

Cone and No cone, T = 900 K 

Mkn 348 

7 

0.4 - 1.0 

36° - 44° 

All 

Mkn 533 

7 

1.5 

30° - 36° 

Cone and No cone, T — 900 K 

Mkn 573 

6 

0.4 - 0.8 

36° - 50° 

Cone and No cone, T = 900 K 

MCG— 5-23-16 

5 

0.7 -0.8 

30° - 36° 

Cone and No cone, T = 900 K 

IC 4329A 

6 

1.4 

10° 

Cone, T = 900 K 


Note. — X 2 /v 1S X 2 P er degree of freedom, where v = n — 2 and n is number of points for 
the SED. 
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NGC106B, Model T=1500K, no cone 





NGC3227, Model T=900K, no cone 


05 



NGC3281, Model T=900K, no cone NGC4151. Model T=900K, cone 




, ^ ,G - * 7 Observed SEDs and best-fitting models for the sample of Seyfert galaxies. The viewing angles (0,.) for the model outputs 
rom top to bottom, and are given in Table 6. For NGC 7469 we show the two best fitting models (cone and coneless). 


increase 
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log yS y (l.6/zm)/i'S^(4.8/xm) 


log i/S l/ (1.6/xm)/i> , S l/ (4.8/zm) 



log i^S l/ (1.6^,m)/uS l/ (3.8/i,m) 


log uS v (1.6fxm)/uS v (3.8fxm) 


Fig. 2. — Non-stellar near-infrared color-color maps. The left panels show the data for the sample of Seyfert galaxies, divided into galaxies 
with viewing angles 0 V < 30° and 0 V > 30°. The right panels show the outputs from Efstathiou et al. (1995) models for the torus+cone 
component model (thick lines), and for the torus alone (coneless) model with two dust sublimation temperatures, T = 900 K (dot-dash lines) 
and T = 1500 K (solid lines). The crosses on the lines indicate the viewing angles (see Table 6) which increase to the left. We also show the 
average position of low redshift PG quasars (see text). 
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this model is /b = 6. 

The viewing angles with 6 V = 0° polar view, and 6 V = 
90° equatorial view (note that in Efstathiou & Rowan- 
Robinson 1995 models the authors measure the angles 
from the equator), UV optical depth, and corresponding 
visual extinctions of the coneless and cone torus model 
outputs are given in Table 6. 

4.2. SED fitting 

The first step was to scale the outputs from the model 
to the observed SEDs in v S v . The SED fitting was per- 
formed by minimizing the x 2 function of the models and 
the observations for the varying viewing angles to the AGN 
and normalizations at different wavelengths, taking into 
account the errors due to the determination of the non- 
stellar fluxes and the photometric uncertainties. Note that 
we have not attempted to fit NGC 2992 or NGC 4968 be- 
cause there are only four points for their SEDs. In Table 7 
we list the best fitting models (at la level confidence) for 
the SEDs where n is the number of points in the SED, 
X 2 /v is x 2 P er degree of freedom, 6 V is the best fitting 
viewing angle (or range of angles) and the corresponding 
models. 

In Figure 1 we present the observed SEDs (logz/S^ in 
units of erg cm -2 s" 1 ) for the sample of Seyferts along with 
the best fitting model scaled to the observed SED. The 
model outputs are plotted for the viewing angles given in 
Table 6 (0 V = 0° is the top line, and 0 V increases to the 
bottom). The abrupt transition seen in the model output 
at viewing angle 0 V = 30° is due to the assumed beaming 
factor. Using a beaming factor varying smoothly with 9 
would avoid the abrupt transition between those viewing 
angles. 

As can be seen from Table 7 and Figure 1, the cone 
model provides good fits to the observed SEDs except 
in the cases of some intermediate SEDs (the SEDs of 
NGC 3227 a M NGC 5506). The reason why the cone 
model fails to reproduce these intermediate SEDs is be- 
cause of the abrupt transition produced for the assumed 
form of the beaming factor. If the beaming factor is de- 
creased, then the cone model produces an acceptable fit 
to the SEDs of these galaxies statistically indistinguish- 
able from the coneless model. In some cases we find from 
our SED fitting that both the cone and coneless models 
produce indistinguishable fits at the la confidence level. 

Despite the fact that we do not find compelling evidence 
to favor either the coneless or the cone model alone, we find 
that the viewing angle to the central source is very well 
constrained by the modelling of the SEDs. We find that 
none of the galaxies (not even the hard X-ray Compton- 
thick galaxies) in our sample show viewing angles greater 
than 64° , with most of the obscured Seyferts displaying 
values of between 30 and 40°. This does not mean that 
there are no Seyfert galaxies with 6 V > 65°. First, the 
values for the fitted viewing angles will depend on the as- 
sumed geometry for the torus model (mainly, the equato- 
rial optical depth and the beaming factor). Second, the 
shapes of the SEDs for viewing angles above 65° become 
almost parallel, with the added complication that the fit- 
ted viewing angle is very sensitive near-infrared non-stellar 
estimates. Finally, we stress that our sample is not com- 
plete or unbiased in the sense that we are most likely not 
including very obscured objects. 


All the Seyfert Is (see also Granato et al. 1997) except 
NGC 3227 show relatively flat SEDs implying viewing an- 
gles in the range 0 V = 0° - 30° - almost a face-on view. 
The coneless torus model does not account for the bright 
1 - 5 ^m emission present in NGC 4151 and IC 4329A. 
Other models also invoke additional near-infrared compo- 
nents in order to explain the SED of NGC 1068. Granato 
et al. (1997) assumed a model similar to the ’anisotropic 
spheres’ of Efstathiou & Rowan- Robinson (1995), whereas 
Pier & Krolik (1993) used a hot component made up of 
clouds inside the inner radius of the torus. 

As expected, the Seyfert 2 galaxies show SEDs with 
different shapes well explained with the varying viewing 
angle to the central source. For example, NGC 5252 and 
NGC 5506 (both classified as Seyfert 1.9) appear to be seen 
at a viewing angle similar to that NGC 3227 (classified as 
a Seyfert 1.5). For the particular model assumed here it 
would be 6 V = 30°. The similarity between the SEDs of 
NGC 5506, NGC 5252 and NGC 3227 is also supported 
by the detection in NGC 5506 and NGC 5252 of a broad 
line component in the infrared hydrogen recombination 
line Pa/? (Rix et al. 1990; Blanco, Ward, & Wright 1990; 
and Ruiz et al. 1994, respectively). Veilleux et al. (1997), 
however, argued against the presence of such broad compo- 
nent in NGC 5506. Both galaxies also show relatively low 
hydrogen column densities derived from hard X-ray obser- 
vations (ATh = 4.3and3.4 x 10 22 cm“ 2 for NGC 5252 and 
NGC 5506 respectively Bassani et al. 1999; equivalent to 
an optical extinction of Ay ~ 22 and 17 mag for the stan- 
dard gas-to-dust conversion factor). NGC 5506, NGC 5252 
and NGC 3227 appear as good examples of why averaging 
properties of type 1 and type 2 objects may result in a loss 
of information. 

The rest of the Seyfert 2s are viewed at angles 6 V > 30° 
(NGC 1068, NGC 3281, NGC 7172, Mkn 348, Mkn 533, 
and MCG— 5-23-16). Of these galaxies Mkn 348, Mkn 533 
and MCG— 5-23-16 show broad components in infrared hy- 
drogen recombination lines (Ruiz et al. 1994 for the first 
two and Veilleux et al. 1997 for MCG-5-23-16). Note 
however, that Veilleux et al. (1997) disputed the finding 
of broad lines in Mkn 533. The optical extinctions derived 
from hard X-ray observations (taken from Bassani et al. 
1999) vary from moderate (MCG-5-23-16, Ay ~ 8 mag), 
intermediate (NGC 7172 Ay ~ 60 mag and NGC 3281 
Ay ~ 400 mag) to Compton-thick galaxies (NGC 1068 
and Mkn 533). The most obscured SEDs in terms of the 
fitted viewing angles are Mkn 1 and Mkn 573. Neither 
galaxy shows broad infrared hydrogen recombination lines 
(Veilleux et al. 1997). 

We find a tendency for the galaxies in our sample to 
have lower fitted extinctions from the SED modelling than 
those derived from hard X-ray hydrogen column densities 
for those galaxies with relatively high values of N\i. The 
particular values of the extinctions derived from fitting the 
SEDs depend on the assumed equatorial UV optical depth. 
As discussed in Efstathiou & Rowan- Robinson (1995) in 
order to reproduce the high extinction needed for Compton 
thick objects, it would be necessary to increase the value 
of t\jv up to 4000 mag. This would cause the models to 
produce a very deep 10 fim absorption feature. The main 
difficulty in tori with optical depths of the order of a few 
thousands, however, is that the SEDs become too narrow 
as the temperature drops more steeply with radius. This 



log f(9.6/im) (erg c 


log f(HX) (erg cm 2 s l ) 


log L(HX) (erg s *) 


Fig. 3 — M-band versus hard X-ray (2 - lOkeV) correlations in fluxes (left panel) and luminosities (right panel) 



log f(HX) (erg cm 2 s *) log L(HX) (erg s l ) 


FlG. 4. — As Figure 3, but with the ISO 9.63 /xm fluxes and luminosities. 
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is evident in the models of Pier & Krolik (1993), and it is 
basically the reason why the ’optically thin dust’ in their 
model for NGC 1068 is cooler than ours. But even if we 
were to assume a much higher equatorial optical depth, 
the infrared SEDs of two Compton thick galaxies in our 
sample (NGC 1068 and Mkn 533) are well reproduced with 
relatively modest extinctions that could not possibly be as 
elevated as those inferred from hard X-ray column densi- 
ties. The relatively low values of the extinction needed 
for NGC 1068 were already noted in previous works (e.g., 
Efstathiou et al. 1995; Granato et al. 1997). We inter- 
pret the differing hard X-ray and infrared extinctions in 
the context of the torus -F cone model. Most of the near- 
infrared emission seems to originate in the extended cone 
configuration which suffers less extinction than radiation 
emitted much closer to the active galactic nucleus (e.g., 
the hard X-rays). 

4.3. Non- stellar color- color diagrams 

A number of works have attempted to use near-infrared 
to mid-infrared colors or flux ratios to constrain the ob- 
scuring torus geometry, to derive viewing angles to the 
torus, or to simply quantify the different observed proper- 
ties between type 1 and type 2 Seyferts. For instance, Pier 
& Krolik (1993) used the L — N colors for Seyfert galaxies 
to constrain the ratio ajh (that is, the ratio of the inner 
radius of the torus to its thickness, and therefore the open- 
ing angle of the torus), but failed to reproduce the L — N 
colors of Seyfert 2s, as well as the flat SEDs of Seyfert Is 
and PG quasars. Murayama et al. (2000) have used Pier 
& Krolik (1993) models and the ratio of the L (3.5 ^m) 
and IRAS 25 //m fluxes to derive torus properties such as 
the opening angle and the critical viewing angle for which 
the broad line region (BLR) is visible. Fadda et al. (1998) 
compared the average of number of near-infrared and mid- 
infrared colors for Seyfert Is and 2s with predictions of a 
wide series of models. These studies, as mentioned above, 
may suffer from unquantified large aperture effects in the 
mid-infrared where the underlying galaxy and/or star for- 
mation may make a significant contribution. 

We explore non-stellar color-color diagrams involv- 
ing observations with A < 5 fim for which we are 
certain about the non-stellar contribution. Due to 
the large apertures employed by the ISO instru- 
ments, there may be a significant stellar contribu- 
tion. In Figure 2 we present non-stellar color-color di- 
agrams: V - M vs. H - M [uS v (3.8 fim)/vS u (4.8 fim) 
vs. vS v ( 1 . 6 /mi) (vS v (4.8 /xm), upper panel] and 

V — M vs. H — V [^5^ (3.8 fim)/i/S l/ {4.8 fim) vs. 
vS v (1.6 fvm)jvS v (3.8 /xm), bottom panel]. The left pan- 
els show the sample of Seyfert galaxies. For the Seyfert Is 
in Table 5 the 3.8 /im fluxes have been interpolated using 
the 3.5 fim and 4.8 fim fluxes. The errors at 1.6 fim account 
for the uncertainty associated with the non-stellar flux de- 
termination, whereas the errors at 3.8 pm and 4.8 pm are 
due to the photometry and sky subtraction uncertainties. 
We have divided the sample of Seyferts into galaxies with 
fitted viewing angles 6 V < 30° and 6 V > 30°. NGC 2992 
and NGC 4968 are plotted with the symbol for unknown 
6 V since they were not fitted in the preceding section. In 
the right panels we show the outputs from the torus + 
cone model, and the coneless torus models with dust sub- 
limation temperatures T = 900, 1500 K. The crosses on 


the lines are the different viewing angles (increasing to 
the left). We also show the average position of nearby 
PG quasars using a near-mid infrared spectral index of 
a — 1.4 ± 0.13 (Neugebauer et al. 1987). 

The effect of adding the cone component is readily seen 
from this figure. The cone model is able to reproduce the 
bluer (that is, flatter SEDs) H - L\ H - M and V - M 
colors of the average PG quasar, and the Seyferts with 
0v < 30°. Note that the three Seyfert Is located in the 
upper right corner (NGC 3227, NGC 4151 and IC 4329A) 
have anomalous M-band fluxes (too low) as seen from the 
SEDs in Figure 1. In addition, the coneless models pro- 
duce much redder colors for large values of the viewing 
angles, but these axe not observed in our sample. 

These two color-color diagrams provide a good means 
to separate Seyfert 2 galaxies with modest obscuration, 
Ay up to ~ 20 mag from hard X-ray observations (e.g., 
NGC 5252, NGC 5506 and MCG -5-23-16, NGC 2992), 
from those with high obscuration. The two galaxies whose 
SEDs were not fitted in the preceding section also fit 
into these two categories. NGC 2992 displays a mod- 
erate absorption from haxd X-ray observations, whereas 
NGC 4968 is a Compton-thick galaxy (Bassani et al. 
1999). There appears to be a transition at around 
log vS„ (3.8 fim)luS v (4.8 fim) ~ —0.15. This may be in- 
terpreted as corresponding to a critical viewing angle for 
a relatively hard edge torus, above which the BLR be- 
comes visible. This critical value is larger than the open- 
ing angle. The separation between the two classes is 
not so well defined if the logi/S„(1.6 fim) /uS^ {4.8 fim) (or 
log vS v (1. 6 fim)/vS l/ (3.8 /im)) flux ratio is used instead, 
probably due to other effects such as foreground extinc- 
tion still affecting the 1.6 fim fluxes. 

5. MID-INFRARED EMISSION OF AGN, AN INDICATOR OF 
THE AGN LUMINOSITY? 

5.1. The mid-infrared versus hard X-ray correlations 

As mentioned in the introduction all torus models pre- 
dict strong emission at mid-infrared wavelengths. Based 
on ISO observations of quasars it has been found that there 
is a good correlation between the infrared (3 - 40 /xm) lu- 
minosities and other properties (blue luminosity, soft X- 
ray luminosities) and this is indicative that an important 
fraction of the the mid-infrared energy is produced by dust 
heated by the AGN (Polletta et al. 2000; Haas et al. 2000). 
To test if this is the case for both Seyfert Is and Seyfert 
2s, it is necessary to have a good indicator of the AGN lu- 
minosity to which we can ratio the mid-infraxed emission. 

The hard X-ray (2 - 10 keV) emission of Seyfert galaxies 
is known to be a good indicator of the intrinsic luminosity 
of the AGN for those cases where it is transmitted through 
the torus, that is, in Compton thin galaxies. Other pro- 
posed indicators of the AGN power of both Seyfert Is and 
Seyfert 2s include the [Om]A5007 luminosity (Mulchaey 
et al. 1994; Heckman 1995) or the non-thermal 1.45 GHz 
radio continuum (Heckman 1995). The [O III] luminosities 
may be affected by extinction, distribution of the gas in 
the parent galaxy, and star formation. Heckman (1995) 
found that the 10.6 fim emission of Seyfert 2s when ra- 
tioed to the non-thermal 1.4 GHz emission is four times 
smaller than that of Seyfert Is, and this would suggest that 
even the mid-infrared emission of Seyfert 2s is anisotropic. 



16 


The non-stellar infrared continuum of Seyfert galaxies 


However, caution is needed when interpreting Heckman’s 
results; the radio emission of Seyfert galaxies may not be 
an isotropic property because of Doppler boosting of the 
core emission in Seyfert Is. Due to the uncertainties as- 
sociated with the use of the radio and [0 in] luminosities 
as indicators of the AGN power, we will test the above 
assumption using hard X-ray measurements. 

We use the M - band (4.8 /xm) measurements presented 
in this work for Seyfert 2 galaxies, and for other Seyferts 
(mainly Seyfert Is compiled by Ward et al. 1987), and the 
ISO 9.63 /xm fluxes from Clavel et al. (2000), and hard 
X-ray (2 - 10 keV) measurements (compiled by Mulchaey 
et al. 1994 and Alonso-Herrero et al. 1997, and refer- 
ences therein) and new measurements from Bassani et al. 
(1999). The hard X-ray fluxes of Seyfert 2s have been 
corrected for the intrinsic absorption. 

In Figures 3 and 4 we plot the correlations in fluxes and 
luminosities (for Ho = 75kms 1 Mpc *) for the 4.8 /xm 
and 9.63 /xm measurements respectively. Both figures pro- 
vide good evidence that the 5 — 10 /xm emission in AGNs 
is powered by the central source. Note that in these two 
diagrams only Seyfert 2 galaxies with relatively low atten- 
uations for which the hard X-rays are transmitted (Comp- 
ton thin, Nn < 10 24 cm" 2 , or Ay <; 500 mag) are plotted. 
A similar correlation has been found by Krabbe, Boker, & 
Maiolino (2000) using N-band (at 10.2 /xm) imaging of a 
small sample of Seyfert galaxies. These two correlations 
between the hard X-ray and the mid-infrared luminosi- 
ties are very much improved even for Seyfert 2 galaxies 
when compared with shorter wavelengths (2.2 and 3.5 /xm, 
Alonso-Herrero et al. 1997). 

We have performed a statistical analysis of the correla- 
tions (both in fluxes and luminosities) using the ASURV sur- 
vival analysis program (Isobe & Feigelson 1990; LaValley, 
Isobe, & Feigelson 1992). The results are tabulated in Ta- 
ble 8. For the Z^band vs. hard X-ray correlations we have 
used the data for Seyfert Is and 2s presented in Alonso- 
Herrero et al. (1997). In this study we demonstrated 
that for some Compton-thin Seyfert 2 galaxies much of 
the L-band (3.5 /xm) emission is attenuated, and this fact 
was used to estimate the extinctions to the AGN. The 
reason for the improved correlations in the mid-infrared 
(4.8 /xm and 9.6 /xm) is the reduced extinction as the wave- 
length increases. As we saw in Section 4.2, in our sample 
there are no Seyfert galaxies with viewing angles 0 V > 64°. 
This means that we should not expect extinctions higher 
than A v = 130 mag or A( 4.8 /xm) ~ 3 mag (using Rieke 
& Lebofeky 1985 extinction law), which is in good agree- 
ment with the observed scatter in the 4.8 /xm vs. hard 
X-ray correlation. 

As discussed in Section 4.1, in the torus -h cone model 
most of the near-infrared emission is produced in the op- 
tically thin cone. This component is located between the 
BLR and the NLR. The conical dust is more visible from 
all viewing angles, and it is probably illuminated by a 
stronger continuum (due to the beaming). This is con- 
firmed by the findings of Bock et al. (1998) that NGC 1068 
shows dust emission from the wall of the cavity, i.e., per- 
pendicular to the axis of the torus, and in an extended 
component (~ 100 pc) coincident with the radio emission 
and ionization cone (Alloin et al. 2000). This naturally 
explains the good correlation between the 4.8 and 9.7 /xm 
and the hard X-ray emissions. Moreover, this may offer 


an explanation for the Compton optically thick Seyferts 
(Nh > 10 24 cm" 2 ) that are bright in the mid-infrared. 
The examples in our sample are NGC 1068, Mkn 533, and 
NGC 4968. This kind of objects are also present in the 
CfA sample (Quillen et al. 2001, in preparation). We con- 
clude that the 5 — 10 /xm emission in Seyfert galaxies is a 
good indicator of the AGN power. 

5.2. Mid-infrared counterparts of highly obscured X-ray 

galaxies 

The use of mid-infrared measurements as indicators of 
the AGN luminosity may have very important implica- 
tions. Recent deep field surveys from Chandra have found 
a number of hard X-ray sources which are not identified 
optically (e.g., Mushotzky et al. 2000). These sources 
are candidates for the highly absorbed objects needed to 
account for the spectral index of the hard X-ray back- 
ground (Setti & Woltjer 1989). Some of the X-ray iden- 
tified sources require high levels of absorption to account 
for the flatness of their X-ray spectrum, suggesting that 
these sources might be opaque even at 10 /xm. A recent 
comparison of faint SCUBA sources with Chandra sources 
found only one source in common (Fabian, Smail, & Iwa- 
sawa 2000). However, in the preceding section we have 
shown that Compton-thick sources can be bright in the 
mid-infrared. If such sources are common then follow up 
deep surveys can search for mid-infrared counterparts to 
highly absorbed X-ray sources. 

For the non-Compton thick sources (Nn < 10 24 cm~ 2 ) 
the 5 /xm luminosity is similar to that at 2 — 10 keV (Fig- 
ure 3). Mushotzsky et al. (2000) found 2000 and 300 
sources per square degree with 2 — lOkeV fluxes greater 
than 2xl0” 15 and 10" 14 erg cm -2 s" 1 respectively. These 
sources would correspond to fluxes of 3 and 17/xJy at rest 
wavelength of 5 /xm. If the rest 10 and 20 /xm luminosities 
are twice that of 5/xm flux then we would expect 13 and 
66/xJy at a rest wavelength of 10 /xm and 26 and 130/xJy 
at a rest wavelength of 20 /xm for these sources. 

The K correction in the mid-infrared is likely to make 
the observed sources more difficult to detect since the spec- 
trum is red out to ~ 30 /xm. The obscured sources are most 
likely to be detected between - 15 - 90 /xm where the rest 
wavelength 5 - 30 /xm peak is redshifted into the observed 
wavelength range, though unobscured sources would suf- 
fer less from the K correction. Deep mid-infrared sur- 
veys planned with SIRTF will go deep enough over large 
enough areas that mid-infrared candidates to absorbed X- 
ray sources should be found. However we expect based on 
ISO surveys that in this wavelength region the number of 
counts from starbursts will dominate by at least an order 
of magnitude the counts from AGNs. If these sources are 
at a redshift of 2.0 then the X-ray luminosity of a source 
with hard X-ray fluxes 10” 14 erg cm" 2 s" 1 is only ~ 10 44 
erg s" 1 (assuming qo = 0.15 and Ho — 75) which implies 
that these sources might just be Seyferts. If they are L* 
galaxies they would be difficult to detect in optical follow- 
up searches (mi > 26). 

6. SUMMARY AND CONCLUSIONS 

We have presented 1 - 5 /xm imaging for a sample of 
Seyfert 2 galaxies. The data have been complemented with 
existing HST/NICMOS (1 - 2.2 /xm) and ISO (6 - 15 /xm) 
observations of the Seyfert 2s, as well as a small compari- 
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Table 8 

Correlation statistics for Seyfert Is and Seyfert 2s. 


Correlation 

quantity 

N 


probability 

L-band vs. HX 

flux 

33 

0.26 

1.4 x 10 -1 

M-band vs. HX 

flux 

24 

0.61 

3.3 x 10" 3 

9.6 fim vs. HX 

flux 

33 

0.66 

2.0 x 10" 4 

L-band vs. HX 

lum 

33 

0.90 

< 7 x 10" 6 

M-band vs. HX 

lum 

24 

0.94 

6.8 x 10~ 6 

9.6 vs. HX 

lum 

33 

0.90 

6.7 x 10~ 7 


Note. — N is the number of objects in the sample, r s is the 
Spearman correlation rank, and prob is the probability that the 
correlation is not present. 


son sample of Seyfert Is. We have performed a very care- 
ful estimate of the non-stellar fluxes in the 1 — 5 fi m range 
to construct non-stellar infrared SEDs. We fit the non- 
stellar infrared (up to 15 fim) SEDs with the Efstathiou 
& Rowan- Robinson (1995) code for tapered disks. We 
consider a simple torus model (coneless) and a compos- 
ite model which includes a torus and an optically thin 
cone component. The latter model was found to fit well 
the SED of NGC 1068 (Efstathiou et al. 1995). From our 
analysis we find that: 

• The non-stellar infrared SEDs for the Seyfert galax- 
ies in our sample display a variety of shapes, which 
can all be well fitted by either cone or coneless mod- 
els. Although the model fits to the SEDs do not 
clearly favor the cone models, as claimed by Efs- 
tathiou et al. (1995), we nevertheless find that the 
viewing angle towards the central source is well con- 
strained even in cases where the fits to the cone and 
coneless models are equally valid. For all the galax- 
ies in our sample we find viewing angles 6 V < 64°, 
with most of the obscured galaxies showing 8 V ~ 40°. 
The fitted viewing angles are however dependent on 
the assumed torus geometry. 

• The non-stellar V — M vs H — V and H — M 
color-color diagrams provide a good means to sep- 
arate those Seyfert 2s with moderate obscurations 
(Ay 20 mag from hard X-ray observations) from 
those with high obscuration. The V — M, H — V 
and H — M colors of Seyfert Is and Seyfert 2s with 
viewing angles of 6 V < 30° are better reproduced 
with the cone model. 

• The extinctions derived from the model fits to the 


SEDs tend to be less than those inferred from mea- 
surements of the hard X-ray attenuations. This can 
be understood if the material responsible for obscu- 
ration of the infrared continuum is in a different lo- 
cation from the gas column causing absorption of the 
X-rays. 

• There is a good correlation between 4.8 fim and ISO 
9.7 fi m and hard X-ray fluxes and luminosities for 
both Seyfert Is and Compton thin (Nu < 10 24 cm -2 ) 
Seyfert 2s. The improved correlations at 4.8 and 
9.7 fim with respect to those at shorter wavelengths 
are explained in terms of the reduced extinction. 
Some Compton thick sources (e.g., NGC 1068 and 
Mkn 533) are bright infrared sources suggesting that 
the component responsible for the bulk of the in- 
frared emission in Seyfert galaxies is more visible 
from all viewing angles than that responsible for the 
hard X-ray emission. We conclude that the mid- 
infrared emission in Seyfert galaxies can be used as 
a measure of the AGN luminosity. This has impor- 
tant implications for future mid-infrared searches of 
visually obscured objects which may make an impor- 
tant contribution to the hard X-ray background. 
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